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Abstract 

We present results suggesting that magnetic monopoles can account for chiral 
symmetry breaking in abelian gauge theory. Full U{1) configurations from a 
lattice simulation are factorized into magnetic monopole and photon contribu- 
tions. The expectation {tpip) is computed using the monopole configurations 
and compared to results for the full U{\) configurations. It is shown that 
excellent agreement between the two values of (^V') is obtained if the effect 
of photons, which "dress" the composite operator iptp, is included. This can 
be estimated independently by measurements of the physical fermion mass in 
the photon background. 
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f/(l) lattice gauge theory provides an ideal laboratory for studying the effects of 
monopoles. In U{1), all non-perturbative effects are known to be caused by monopoles 
Further, there is an efficient procedure for identifying the magnetic current of monopoles in 
the link angle configurations generated in a U{1) lattice gauge theory simulation [H. Pre- 
viously, two of the present authors have studied the heavy quark potential in f/(l) lattice 
gauge theory, and shown that in the confined phase, the value of the string tension calculated 
from monopoles agrees quantitatively with the full U{1) string tension calculated directly 
from link angles 0. 

In this paper, we turn our attention to the quenched chiral condensate in U{1) lattice 
gauge theory, (ipip). We expect (and find) that this has a non-zero value in the confined 
phase, and vanishes in the deconfined phase. Our interest is in how the value of (ipip) cal- 
culated from the link angle configurations compares with the same quantity calculated from 
monopoles. Ascribing non-perturbative effects to monopoles would say that in the confined 
phase, a non-vanishing condensate will be induced by monopoles. Nevertheless, since a ma- 
trix element of a product of field operators is involved, this value should be renormalized 
by short distance perturbative effects i.e. photons, so to make a quantitative comparison, 
we also need to calculate this renormalization. Our method for doing so, described in more 
detail below, is to compute the ratio of renormalized to bare masses for a charge moving 
solely in the field of photons. Our principle result is that applying this renormalization 
factor to the monopole value of (ipip) yields the full U{1) answer, to within statistical errors. 

Let us now briefly describe the factorization of the link variables into photon and 
monopole parts. We start by resolving the U{1) plaquette angles (j)^u into fluctuating and 
monopole parts, 

(Pi^u = + 27rm^^, 

where G (— vr, +7r]. From the monopole term m^j,, we can define two currents, one 
electric, the other magnetic. The electric current is defined by = V~m^jy. The magnetic 
current of monopoles is defined by = V^*m^^, where *m^i, is the dual of m^^. The 
currents i^{m^) reside on direct(dual) lattices. 

We will compute the chiral condensate using the Banks-Casher formula |^ (see below) 
which involves the density of eigenvalues of the Dirac operator. Since the Dirac particle is 
electrically charged, the effect of monopoles on it must be represented by an electric vector 
potential. This four- vector potential derives from the electric current j^, and is analogous to 
the familiar three- vector potential A used for the case of a static monopole in the continuum. 
We have 

A^"-{x)^gY.v{x-y)3,{y\ (1) 
y 

where g = 27r/e is the magnetic unit of charge, and v the lattice Coulomb propagator in 
Feynman gauge, satisfying V^V~u(x — y) = —Sx,y By defining UJ^""'{x) = exp{ieA'^""-{x)), 
we define the photon link by demanding that {f/^} be factorized on each link: 

U,ix) ^ Ur{x)Uf°\x). (2) 

We can now examine the chiral condensate {ipip) separately on configurations {t/^}'""" 
and {f/^}'''^"*, as well as the full {U^}. In the quenched approximation {ipip) can be defined 
in the chiral limit via the eigenvalues of the Dirac operator 0]: 
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m = T7P(A = 0), 



(3) 



where V is the lattice volume and p{X)dX is the number of eigenvalues satisfying p\n) = 
iXn\n) in the interval (A, A + dX), with p[Uf^] the staggered lattice fermion kinetic operator 

PM = I Y.i-^r^'"^^''-' (t^.(^) W - u;{y)Sy,.-f:) . (4) 

For finite V, p(0) vanishes; therefore an extrapolation to the A — > limit must be made 
from the spectrum density at non-zero A. In practice p(A) can be estimated by measuring 
the lowest 0(25) eigenvalues of p per configuration using the Lanczos algorithm, and then 
using a binning procedure 0. 

Calculations of the chiral condensate were performed with U{1) configurations generated 
using the standard Wilson action on a 12^ lattice. Three different values of the inverse 
coupling P were chosen: P = 1.005, 1.010, and 1.020. For each value of P, 5000 lattice 
updates were performed before beginning measurements to allow for equilibration. The next 
4000 configurations were used for measurements skipping every 20 lattice updates to reduce 
correlations. For values /? < 1.010, our results exhibited the breaking of chiral symmetry 
(i. e. {ijjijj) 7^ 0). These two values of jS are known to correspond to the confinement phase 
as well 0. For /5 = 1.020 the configurations were found to be in the chirally symmetric 
phase, which is coincident with the Coulomb phase. 

Our results for the chiral condensate calculations for the chirally broken phase are pre- 
sented in Figs. (1) and (2). We show both the full U{1) calculations and the results from the 
monopole gauge field configurations obtained using Eq. (1). Using a linear fit and Eq. (3), 
the value of the chiral condensate was extracted from the 10 lowest values of p(A). The 
values for (ipip) are shown in Table I. 

In Fig. (3) the results for the chiral condensate from full U{1) fields and monopoles are 
shown for /3 = 1.020. For comparison, the results from the /? = 1.010 monopole configura- 
tions are also included. Linear fits to p{X) yield a very small, but finite intercept. However, 
the values are found to be about a factor of 100 times smaller than the results in the broken 
phase. Thus, we are confident of being in the chirally symmetric phase for /5 = 1.020. For 
completeness, the P = 1.020 results are also included in Table I. 

In Fig. (4) we show the eigenvalue spectrum calculated using the {U^}^^"'^ background 
from both broken and symmetric phases. There is no signal of chiral symmetry breaking 
(note that we choose anti-periodic boundary conditions for the fermions in the temporal 
direction to avoid zero modes from near-plane wave solutions). This supports our conjecture 
that chiral symmetry breaking in this model can be ascribed entirely to monopoles. 

Is it possible to account for the mismatch between {xjjxjj)^^^^ and {ipip)"^°"7 Since ipip is a 
composite field operator, we expect it to be modified by quantum corrections independently 
of whether it acquires a vacuum expectation value. In perturbation theory, UV fluctuations 
in general result in the requirement for all field operators to be renormalized, and for com- 
posite operators to have an additional renormalization. This consideration leads us to the 
following hypothesis: the mismatch between the two condensate measurements is due to the 
rescaling of the iptp operator by the fluctuations of the gauge fields contained in the {f/^}^^°* 
configurations; more concisely 
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= ((^7/^)^'^°*)'"°", (5) 

where [ipipy^"^ is the dressed operator which takes into account photon-hke fluctuations, 
and the expectation value is taken in the monopole-only backgrounds. If the hypothesis is 
correct, then we can write 

(^^)^'°* = Z{^P^pr^'^, (6) 

ie. the tree-level operator is multiplicatively enhanced by the photon fluctuations, the con- 
densate data suggesting that the factor Z ~ 1.5. 

In support, we now describe an alternative and independent determination of Z, obtained 
by measuring the physical fermion mass in the {t/^}^'^"* background. On the assumption 
that the fluctuations in {f/^}^'^°* are approximately Gaussian, then the dressing of the ipil) 
operator is given by a set of Feynman diagrams. Let denote the complete set of IPI 
diagrams describing corrections to the fermion two-point function. Then for the dressed 
fermion propagator we have 

tp + mo — ^{p) + rriR 

where mo, denote bare and physical fermion masses respectively, Z2 is a wavefunction 
rescaling which by analogy with perturbative QED we expect to be both gauge and cutoff- 
dependent, and the ~ symbol shows that the second equality holds only in the neighborhood 
of some subtraction point ip = fi. Equation (|^) may be rearranged to read 



with 



Now, using the identity 



S(^) = (1 - Z^'Z,n)mo + - Z,') (8) 



niR = Zm{a, ij)mo. (9) 



d 1 _ -1 

dmQ ip + niQ {ip + ■mo){ip + tuq) ' 

we see that the operation of differentiating with respect to mo is equivalent to a zero mo- 
mentum insertion of a ipip operator in a Feynman diagram. Therefore —dTi/diriQ is the set 
of IPI diagrams, having one external ip and one external ip^ which describe corrections to 
■xpil), ie. 

(^^)?"^°* = zJi- {^/j^/jy^"' = zi^pi,)'-"-', (11) 



le. 
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where Z{a,fi) is the rescahng factor we seek, and the factor Z2 is included because in the 
diagrammatic approach matrix elements are evaluated with dressed fermion propagators on 
the external legs. Combining with (]T2|), we find 

Z = Zjn = • (13) 

mo 

Hence a measurement of mn in the {f/^}^'^"* background as a function of mo gives an inde- 
pendent estimate of Z. 

This argument is essentially the same as that used to establish Z^^ = mo/m/j in standard 
renormalized perturbation theory, where = Z^^ltpip) is the renormalized operator 

whose Green functions with other renormalized fields and operators are finite [^. Here we 
make no attempt to define a renormalized operator, since bare quantities are evaluated in 
a lattice simulation, but rather use the same formalism to quantify the effects of operator 
enhancement by quantum fluctuations. Note that we have not specified the subtraction point 
defining Z very precisely, and have assumed that Z is independent of mo. From experience 
in perturbation theory, we expect that numerically the most significant contribution to Z, 
of O(lna), comes from short wavelength fluctuations, and that Z is relatively insensitive to 
the details of the subtraction. 

To measure the physical fermion mass in the background, we performed cal- 

culations of the fermion propagator starting with the same configurations used to compute 
p(A). Since the fermion propagator is not gauge invariant, it is first necessary to fix a gauge. 
Although a gauge transformation of {f/^}^'^°* is not strictly a symmetry of the full theory, 
this procedure is justified since ipip, and by hypothesis mn, are gauge invariant. In this work 
we have used Landau gauge [|| and extracted the lattice fermion timeslice propagator 

/ \ phot 

Cf{x,) = Re Yl {mU,]+mo)oi) (14) 

using a conjugate gradient routine, for bare mass values moa = 0.1, 0.09, 0.04. The 
restriction to spatial sites an even number of lattice spacings from the origin in each direction 
improves the signal Note that we have not attempted to fix the residual gauge freedom 
as in 1^, but have instead relied on the fluctuations in {f/^}^'^°* being intrinsically small 
(see below). The physical mass m^j can now be estimated by fitting C/ to the following 
functional form: 

Cf{x4) = A (expi-mnXi) + {-ly* exp(-mjf (L - 0:4))) , (15) 

where L is the lattice size in the time direction. We used 100 configurations from each of 
the three f3 values previously studied. In order to take into account correlations between 
renormalized mass estimates due to using the same configurations for different bare masses 
we applied a bootstrap fitting routine. 

To determine Z we plotted the mass obtained from Cj against the input bare mass, and 
fitted the results to a linear form: 

'^i?("^o) = Z mo + b. (16) 

Fig. (5) shows a graph of m/j(mo) from our simulation for (3 = 1.010. The graphs for the 
other two values of (3 are very similar. The results for Z and b from the linear fits for each 
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value of (3 are shown in Table II. From the table it can be observed that the value of Z is 
a smooth function through the phase transition. Unexpected in our results were the small 
non-zero values of the intercept. This will be discussed briefly below. 

Finally, we compare our value of Z calculated using with the value observed 

in the chiral condensate calculations. We define the renormalization factor from the chiral 
condensate to be 



and present the comparison in Table III. The two independent determinations of the renor- 
malization factor show excellent agreement. This answers the question asked earlier: Is it 
possible to account for the mismatch between (■?/'■?/') and {ijjip)'^°'^7 Our results demon- 
strate that the mismatch is simply operator renormalization due to the photons. 

As yet we have no satisfactory explanation for the non-zero intercepts in Table II, im- 
plying a small breaking of chiral symmetry in the photon-only configurations. It could be 
either a finite volume effect, or perhaps a residual gauge freedom associated with Gribov 
copies and the consequent non-uniqueness of the Landau gauge, or perhaps due to spatially 
constant gauge field modes. The only effect that we have investigated quantitatively is that 
of constant modes. We repeated our previous calculations at /3 = 1.005, this time perform- 
ing a global gauge transformation to remove constant modes from the gauge fields. This 
was done in the following way: The full gauge fields were shifted by a direction dependent 
constant, 

A'^ix) = A^{x) + c^, (18) 

where is chosen so J^x^'^ii^) = 0. The new gauge fields A'^ were used to recompute Z 
and the chiral condensates. In repeating the calculations with this gauge fixing we observed 
no measurable difference in the results. 

Although we are not able to explain the small non-zero intercepts, the overall agreement 
of the slope with the value of Z required to explain the condensate results is extremely satis- 
fying, and supports our hypothesis on the role of the residual photon fluctuations enhancing 
the tpip operator. 

An obvious direction in which to extend this analysis is the exploration of chiral symmetry 
breaking in non-abelian gauge theories following the identification of monopole networks 
after abelian projection. Studies of this kind have already appeared |]TU[; it is interesting to 
note that in the maximal abelian gauge used, there is a similar mismatch in the measured 
values of (ipip) between full and monopole-only configurations. It would be interesting to 
check whether this effect could be accounted for by the operator renormalization described 
here. 

Finally note that even though this argument assumed a split of the background gauge 
field into Gaussian fluctuations, described by Feynman diagrams, and non-perturbative 
monopole-only configurations, the idea of classifying effects into those which rescale a local 
operator and those which lead to a non-vanishing expectation value for the operator may 
be generalized. For instance, it may in principle be possible to extend the analysis by fac- 
toring configurations in a scale dependent fashion, including only large monopole loops in 
the monopole part, and seeing if small monopole loops simply renormalize iptp by comparing 
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the Z factors measured in two different ways. This may prove to be an effective probe of 
the scale at which non-perturbative effects manifest themselves. The scale dependence of 
monopole contributions to the string tension was studied in |]lT . 
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TABLES 



TABLE L The values of {iptp) = yp{0) obtained from linear fits to p{\). 



p 






L005 


0.165(3) 


0.105(1) 


LOlO 


0.138(2) 


0.089(2) 


L020 


0.002(1) 


0.0009(6) 


TABLE IL Results of the fits m/j (mo) (cf. (|T|)) 


P 


Z 


b 


LOOS 


1.570(4) 


0.0124(3) 


LOlO 


1.561(5) 


0.0126(3) 


L020 


1.524(4) 


0.0115(3) 



TABLE in. Comparison of our two determinations of the renormalization factor Z; The column 
labeled Zm comes from the photon configurations via Eq. (|16| ) , and the column labeled Z^^ comes 
from the chiral condensates via Eq. (^) 



p 






1.005 


1.570(4) 


1.570(30) 


1.010 


1.561(5) 


1.550(40) 
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FIG. 1. The spectral density function for /3 = 1.005 from full U{1) fields (solid squares) and 
monopoles (open squares). 
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FIG. 2. The spectral density function for /3 = 1.010 from full U{1) fields (solid squares) and 
monopoles (open squares). 
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FIG. 3. The spectral density function for /3 = 1.020 from full U{1) fields (solid squares) and 
monopoles (open squares). For comparison the results from monopoles for (3 = 1.010 (open circles) 
are also plotted. 
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FIG. 4. The spectral density function from the photon background {11^}^^°^ for (3 = 1.005 
(sohd squares) and /3 = 1.020 (open squares). 
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